Multiple-transit paths and density correlation functions in PASEP 
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We consider the partially asymmetric simple exclusion process (PASEP) when its steady-state 
probability distribution function can be written in terms of a linear superposition of product mea- 
sures with a finite number of shocks. In this case the PASEP can be mapped into an equilibrium 
walk model, defined on a diagonally rotated square lattice, in which each path of the walk model 
has several transits with the horizontal axis. We particularly show that the multiple-point density 
correlation function in the PASEP is related to the probability that a path has multiple contacts 
with the horizontal axis from the above or below. 

PACS numbers: 02.50.Ey, 05.20.-y, 05.70.Fh, 05.70.Ln 



I. INTRODUCTION 



The Partially Asymmetric Simple Exclusion Process, 
or PASEP for short, usually refers to a system of iden- 
tical classical particles moving on a discrete lattice. In 
this system the particles are injected and extracted from 
both ends of the lattice while hopping to the left and 
right in the bulk of the lattice with unequal rates x and 
1 respectively. The particles perform continuous time 
random walks on the lattice. The exclusion rule forbids 
two particles to occupy a single lattice site at the same 
time 

The steady-state of the PASEP has been widely studied 
in the related literatures. The exact results have been ob- 
tained using a matrix product approach for the PASEP 
with open boundaries [j, 0) and also on a ring in the pres- 
ence of a second class particle which is sometimes called 
a defect particle [5*]; however, not much is known about 
the dynamics of shocks which are defined as discontinu- 
ities in the density profile of the particles in the PASEP. 
In the authors have studied the time evolution of 

a product measure, as an initial probability distribution, 
with a finite number of shocks. It has been shown that 
at any time t > Q the probability distribution function 
can be written as a linear combination of shock measures 
with an structure similar to that of the initial configu- 
ration. On the other hand, the shock positions perform 
continuous time random walks interacting by the exclu- 
sion rule. The Monte-Carlo simulations for a single shock 
confirm these results [9, 10] . 

It can be simply understood that if the dynamics of the 
shock position in a one-dimensional driven-diffusive sys- 
tem is similar to that of a simple random walker which re- 
flects from the boundaries of the system, then the steady- 
state of the system can be written in terms of a linear 
superposition of these product measures. On the other 
hand, it has been shown that the steady-state of such 
system can be written in terms of a matrix element of a 
product of non-commuting operators which satisfy an al- 



* Electronic address: ' farhad@ipm.ir| 



gebra. According to this approach, known as the matrix 
product method, one associates an operator to each state 
of a lattice site. The steady-state of the system in this 
case is written in terms of a matrix element (for the open 
boundaries conditions) or a trace (for the periodic bound- 
ary conditions) of a product of such non-commuting op- 
erators which might have finite- or infinite-dimensional 
matrix representations . A review of this method can be 
found in [ll|. 

It is known that if the steady-state of a one-dimensional 
driven-diffusive system with open boundaries can be 
written in terms of a linear superposition of product mea- 
sures containing a single shock, then the steady-state 
of the system can be equivalently expressed by using 
the matrix product method with two-dimensional ma- 
trix representations. It has also been found that the 
shock characteristics appear in the matrix elements of 
the matrix representations i.e. the matrix elements are 
functions of the shock hoping rates and the densities of 
the particles on the left and the right hand sides of the 
shock position Q. As the authors in _8] have proposed 
that this idea can be easily generalized to the PASEP, 
where an A'^-dimensional matrix representation for the 
quadratic algebra of this system exists, which is surpris- 
ingly valid under the same constraints under which the 
steady-state of the system can be written as a superposi- 
tion of product measures with A^ — 1 shocks [6] . The main 
difference between this matrix representation with those 
proposed before is that it is clearly written in terms of 
the hoping rates of the shock fronts and also the densities 
of the particles separated by the them. 
The steady-state of the PASEP has been recently stud- 
ied from different viewpoints. In [l^ a combinatorial 
derivation and interpretation of the weights, associated 
with the stationary distribution of the PASEP, has been 
introduced. By using this approach the authors have 
found explicit expressions for the stationary distribution 
and normalization using both recurrences and path mod- 
els. In a recent work the authors in [isj have used 
a continued fraction representation of the lattice path 
generating function to discuss the steady-state proper- 
ties of the PASEP. In this case an infinite-dimensional 
tridiagonal matrix representation of the quadratic alge- 
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bra of the PASEP has been used to calculate the grand- 
canonical normalization of the stationary distribution of 
the system. It has been shown that by truncating these 
continued fractions one can recover the results obtained 
from a finite-dimensional matrix representations of the 
quadratic algebra which are valid only jilong special lines 
in the phase diagram of this system 14] . 
Apart from these approaches one can try to connect the 
matrix representation of the algebra of a one-dimensional 
reaction-diffusive system directly to the transfer matrix 
of an equilibrium two-dimensional walk model. In this 
way the one can actually show that the normalization 
coefficient of the stationary distribution function of the 
non-equilibrium system coincides with the partition func- 
tion of an equilibrium walk model [isl-flTj . 
In this paper we consider the PASEP with open bound- 
aries where the steady-state of the system can be written 
in terms of a linear superposition of product measures 
with a finite number of shocks. This means that the 
quadratic algebra of the system has a finite-dimensional 
representation under some constraints on the reaction 
rates. It is known that the iV-dimensional representa- 
tions of the stationary algebra describe the stationary 
linear combination of shock measures with A'^ — 1 consec- 
utive shocks Q- Considering an A''-dimensional matrix 
representation of the quadratic algebra we show that the 
partition function of the PASEP is equal to that of a two- 
dimensional walk model on a diagonally rotated square 
lattice. The partition function of the walk model is the 
sum of the weighted paths which start and end on the 
horizontal axis and cross it A^ — 1 times. The multiple- 
point density correlation functions of the PASEP can be 
written in terms of the equilibrium characteristics of the 
walk model. By connecting the finite-dimensional ma- 
trix representation of the PASEP to the transfer matrix 



of the weighted walk model, we can easily calculate the 
multiple-point density correlation function of the PASEP 
in terms of physical quantities defined for the equilibrium 
walk model. 

In order to write a self-explanatory paper we have orga- 
nized it as follows: first we will start with the definition 
of the PASEP with open boundaries and explain how the 
stationary distribution of this system can be written us- 
ing a matrix product method. We will then introduce the 
multiple-transit walk model and find its partition func- 
tion using a transfer matrix method. The multiple-point 
density correlation functions of the PASEP are explicitly 
calculated at the end of the manuscript. 



II. THE PASEP 

We consider a discrete lattice of length L. Each lattice 
site can be occupied by one particle or is empty. Occu- 
pying a single lattice site with two particles is forbidden. 
In the bulk of the system each particle hops to the left 
(right) neighboring site with the rate a; (1) provided that 
it is not already occupied. The particles can enter into 
the system from the leftmost (rightmost) lattice site with 
the rate (1 — x)a ((1 — x)5). The particles can also leave 
the system from the leftmost (rightmost) lattice site with 
the rate (1 — x)^ ((1 — a;)/?). The above-mentioned pro- 
cess, called the PASEP, has been studied widely in re- 
lated literatures. 

Defining a shock as a sharp discontinuity in the density 
profile of the particles on the lattice, it has been shown 
that a product shock measure with A^ — 1 consecutive 
shocks defined as follows 
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can evolve in time according to (A'' — l)-particle dynamics 
provided that the reaction rates lie on a specific surface 
in the parameters space defined by a, /3, 7, S and x. 
The dynamics of each shock position in this case will 
be similar to that of a simple random walker moving 
on a discrete lattice which reflects from the boundaries 
with special rates. In Fig. [T]we have sketched an initial 
distribution of particles on the lattice with A^ — 1 shocks 
locating at different lattice sites. We have assumed that 
the I'th shock hops to the left (right) with the rate Si^i 
(Sr^i) for 1 = 1, ••• ,A^ — 1. In this case the densities and 
the shock hopping rates satisfy the following relations 
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for i = I,-- - ,Ar - 1. (3) 



Using the matrix product method the stationary proba- 
bility distribution vector \P*) for the PASEP on a lattice 
of length L can be written as 



\n = ^{{w\(^\) \v)) 



(4) 
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in which we have defined C := Dq + Di. The operator 
Dq stands for the presence of a vacancy and the operator 
Di stands for the presence of a particle at each lattice 
site. These operators besides the vectors {{W\ and \V)) 
satisfy the following quadratic algebra 



FIG. 1: Sketch of the particle distribution containing N — 1 
shocks. 



DiDo - xDoDi ^ ^{1 ~ x){Do 
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in which Z is a normalization factor and can be obtained 
easily using the normalization condition 
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in which ^ = P7v(l — Pn)- It can be easily verified that 
the following iV-dimensional matrices 
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besides the vectors 
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satisfy ([6]) provided that we have 
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and 
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The boundary rates should also satisfy the following con- 
straint 
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in which we have defined 
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Using the boundary conditions one finds 
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It is interesting to study the properties of the matrix C 
in this representation which is given by 
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It can be seen the ith diagonal element of the matrix C 



is equal to 



pw(i-pw) _ "rr-'^"^ - 
Pi(i-Pi) 



Y\j=i jfr- Considering the case 



< Pi < P2 < ■ ■ ■ < Pn ^ ^, the reader can convince 
himself that only two different cases (phases) exist in the 
steady-state. In the first phase we have 



Pw(l— Pw) 

Pl(l-Pl) 



> 1 



and in the second phase we have ^"(^ P") < i. On the 

^ Pi(i-Pi) 

coexistence line one has ^"^^"'^"^ = i . it is known that if 

pi(i-pi) 

one starts with the product shock measure H]), the shock 
positions move with constant speed until the two shocks 
meet and then coalesce into a single shock. In a long-time 
limit onl y o ne shock (the leftmost, which is the fastest) 
survives [1^. The particle densities on the left- and the 
right-hand sides of this shock are equal to pi and pN re- 
spectively. The quantity gives the ratio of the 
hopping rate of the shock position to the right and left. 
In the first (second) phase the shock position has more 
of a tendency to move to the right (left), and therefore 
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FIG. 2; A typical multiple-transit path containing 5 Dyck 
paths and 4 transits. 



the bulk density will be equal to pi {pn)- 
In the next section we will introduce an equilibrium 
multiple-transit walk model and show how it can be con- 
nected to the PASEP in the above-mentioned case. 



III. THE MULTIPLE-TRANSIT WALK MODEL 

We consider a two-dimensional walk model on a 
diagonally-rotated square lattice. A random walker 
moves along a path which starts at (0, 0) and ends at 
(2L, 0) (see Fig. (2). The random walker starts from (0, 0) 
and is allowed to move one step in the North-East (NE) 
or in the South-East (SE) direction. It can move at most 
two consecutive steps to the NE or to the SE only when it 
crosses the horizontal axis. All other movements are pro- 
hibited. We study the case where the random walker can 
cross the horizontal axis at most N — 1 times. In this case 
the paths taken by the random walker might cross the 
horizontal axis m times where m = 0,-- - ,A'^ — 1; there- 
fore, the path contains m-\-l Dyck paths. We label each 
Dyck path with / where Z = 1, • • ■ ,N. The lih Dyck path 
lies above (below) the horizontal axis if I is odd (even). 
We assign an statistical weight to each path (which might 
contain different combinations of Dyck paths with differ- 
ent lengths) taken by the random walker. The partition 
function of this multiple-transit walk model is then the 
sum of these statistical weights. The statistical weight 
of a path is equal to the multiplication of the weights of 
the individual Dyck paths times the weight of the first 
step times the weight of the last step. In order to find 
the statistical weight of a path of length 2L with m tran- 
sits, the following steps should be taken. The statistical 
weight associated with the lih. Dyck path is defined as 
follows: if I is an even number, we assign the fugacity 
^ to each contact with the horizontal axis except to the 
last upward step, and if I is an odd number, we assign 
the fugacity ^ to each contact with the horizontal axis, 
except the first upward step. If the path starts with the 
1st Dyck path, we assign the fugacity A to the contact of 
the first step with the horizontal axis (equivalently the 
point (0,0)); otherwise, if the path starts with the Zth 
Dyck path [l ^ 1), two different scenarios might happen: 
it is either the first step of the Zth Dyck path (which 
means that it comes from a transit) or the second step 
and above. For the first case we assign a fugacity -i- 



to the contact of the first step with the horizontal axis 
(equivalently the point (0,0)) or J- for the latter case. 

We always assign a fugacity 4- to the contact of the last 
step with the horizontal axis if the path ends to the Zth 
Dyck path. 

Defining the weight of a path as described above one finds 
(after some straightforward but lengthly calculations) the 
partition function of the walk model as follows: 
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(14) 

Up to this point we have not considered any constraints 
on the fugacities. First of all one can easily see that the 
fugacities associated with the boundaries i.e. Si's and z-'s 
do not play any major role in the thermodynamic limit 
(L — > oo) of the partition function ()14|) : however, one can 
consider the cases in which either zi or is the smallest 
bulk fugacity and all other fugacities lie between them. 
This results in two different behaviors for the partition 
function of the walk model in the thermodynamic limit. 
In the first (second) case where zi (z^) is the smallest 
bulk fugacity, the majority of paths contain the Dyck 
path of the first (last) type I = 1 {I = N). This defines 
two different phases which will be discussed later. In 
what follows we will show how the partition function (jl4p 
can be obtained using the transfer matrix method. 



IV. THE TRANSFER MATRIX METHOD 

In this section we will show that the partition function 
of the walk model (fT4)) can be calculated using a transfer 
matrix formalism. By defining a two-step transfer ma- 
trix T = T°T'', in which T° and T" are associated with 
odd and even steps taken by the the random walker re- 
spectively, we will show that the partition function of the 
model (fH|) can be written as 



Zl,n = 



(15) 



One can easily find the transfer matrix T and the vec- 
tors {L\ and \R) as follows: we first define a complete 
base vector associated with the vertices of every path. 
The base vector |0') is assigned to a contact point with 
zero height. For the lih Dyck path, the contact point 
will be from the below (above) the horizontal axis, if I 
is even (odd). On the other hand, the base vector |l') 
is assigned to a point with a unit height either above or 
below the horizontal axis. For the lt\i Dyck path if I is 
even (odd), then the point lies below (above) the hori- 
zontal axis. An example is given in Fig. |31 These base 
vectors are complete in the following sense 
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FIG. 3: Assigning the base vectors to the vertices. 



the base vectors 



= |o') = 



for l = l,--- ,N 



(17) 



in which I is an x identity matrix. For instance in vifhich only the Zth element is non-zero and equal to 
one can choose the following matrix representation for one. If we write the partition function (I14p as 



N N 



Zl^n = 51 ■ ■ • H (i|0'0(0'MT'°|l'')(l''|rlO'^)(0'^| • • • |0'^+0(0'^+^ \R) 



h = l Il + i = 1 



(18) 
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then according to the definition of our walk model the 
only non-zero transitions are given by 



(o'|r°|i'')-<^u' + <5i,i'-i, 



(L|0') 
(0'|i?) 



(19) 



in which 6ij is the usual Kronecker delta. These relations 
imply that we should have the following matrices 
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and vectors 
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(20) 



(21) 



It is now easy to verify that using these matrix represen- 
tation, besides the transfer matrix representation of the 



partition function p5|) . one obtains p4|) . 
An interesting physical quantity which can be defined is 
the probability that the Zth Dyck path has a contact at 
site 2i with the horizontal axis given by 



{L\T'PlT^-'\R) 



{L\Tf\R) ' 

in which the contact operator is defined as follows 



(22) 



Pi = |l')2^-l(0'|2^. 



Normalization requires us to have 



N 



(23) 



(24) 



1=1 



One can easily generalize this idea to find the probability 
that m Dyck paths hjh,--- ,lm {h < ■■■ < In) have 
contacts with the horizontal axis at sites 2ii, 2^2, • • • , 2im 
(«i < • • • < «m). The associated contact operator in this 
case can be written as 
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(25) 
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Now the appropriate probability is given by 



in which 



N 



{L\T^\R) 



(26) 



(27) 



We have assumed that either zi or zn is the smaUest bulk 
fugacity of the walk model. The thermodynamic behav- 
ior of (f26| will be then interesting under these conditions. 
After some calculation one finds in large L limit 



(l-{- 



(28) 
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1 (ZN Ul ' 



n 



k=i ^h,N for > 1. 



It can be seen that these probabilities have exponential 
behaviors with a correlation length which is determined 
by zi and zjv- As explained above, for < 1 {-^ > 1) 
it is more probable that the paths contain the first I = 1 
(last I = N) Dyck path. On the coexistence line where 
zi = zn one finds that the paths might contain both the 
first I = 1 and the last I — N Dyck paths of any length 
not larger than 2L. On this line and in the large L limit 
the probability (1261) is given by the following expression 
which is obviously linear in terms of the position of con- 
tacts 



fc-i 



.N 



(29) 



for fc = 1, • • • , m-|- 1 assuming that io — and im+i — L. 
This expression indicates that on the coexistence line all 
the probabilities are zero unless the paths only contain 
two types of Dyck paths I = 1 and I — N and that the 
first fc — 1 contacts belong to the first Dyck path and the 
last m — k + 1 contacts belong to the last Dyck path. 
In the next section we will show how the multiple-transit 
random walk model is related to the PASEP and that 
how the multiple-point density correlation functions can 
be written in terms of the physical quantities of the walk 
model (USD and 



V. DENSITY CORRELATION FUNCTIONS IN 
THE PASEP 

The reader can easily see that the transfer matrix T in 
has the following N x N matrix representation 
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Comparing this with (1131) one can recognize that they are 
the same matrices provided that 



1 



pjv(l - Pn) 
P»(l - Pi) 



(31) 



for i = 1, • • ■ , iV — 1 and that zj^ — 1. On the other hand 
the partition function of the PASEP ([S]) will be equal 
to that of the multiple-transit walk model (fT4|) provided 
that Vi = Y s-iid = i for i = 1, • • • , iV. In this case 
one has C = T, \V)) = \R) and {{W\ = 
This mapping allows us to write the physical quantities 
of the two systems. For instance the mean density of the 
particles at a given site i in the PASEP {pijLi written in 
terms of the matrix product formalism as 



{Pi)L = 

can be rewritten as 



{{W\C'-^DiC^-'\V)) 

WWW)) 



N 

{Pt)L ^"^Pli 

1=1 



(32) 



(33) 



by noticing that Di = J2i=i Pi'^-Pi- Generalizing this to 
the m-point density correlation function for the PASEP 
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one finds 

N 

{P^,■■■P^JL^ J2 Ph---Pl^{Pnt'''t)L- (34) 

Considering tlie tliermodynamic behavior of (j26p given 
by ((28| and (p9|) one can calculate the thermodynamic 
behavior of the multiple-point density correlation func- 
tion (|34|) for the PASEP. In the same way one can also 
calculate the connected multiple-point density correla- 
tion functions. It can be seen that for the case of two 
point correlation functions the results of [19] can be re- 
covered. For instance on the coexistence line one finds 

rn+l 

{pL., ■ --Pl.Jl - - ^fc-i)pi'~Viv"-'+' (35) 

k=l 

where we have defined Xk — ^ and this is exactly the 
same expression as calculated in po| . 

VI. CONCLUDING REMARKS 

Recently, connections between some of the one- 
dimensional driven-diffusive systems with the equilib- 
rium walk models have been under intensive investiga- 
tions. In this direction it is interesting to show how the 
physical quantities in equilibrium systems are related to 
those in non-equilibrium systems. 



In this paper we have shown that the multiple-point den- 
sity correlation functions in the PASEP are explicitly 
related to the probabilities of having multiple contacts 
with the horizontal axis of a path in an equilibrium two- 
dimensional walk model. In order to show this relation we 
have adopted a matrix product approach while using an 
special finite-dimensional matrix representation for the 
quadratic algebra of the system which explicitly reveals 
the shock characteristics of the PASEP. This matrix rep- 
resentation helps us connect the partition function of the 
PASEP to that of a multiple-transit walk model obtained 
using the transfer matrix method. 

There are a couple of examples in related literature which 
show that the stationary distribution of the asymmetric 
simple exclusion process (ASEP) with open boundaries 
can be related to that of some path models. It has been 
shown that the partition function of each path model, 
obtained using the transfer matrix method, is equal (or 
at least proportional) to that of the ASEP, calculated us- 
ing the matrix product method. 

As we have seen in this paper, some of the physical quan- 
tities defined in both models are related to each other us- 
ing these methods. However, it is still not clear why and 
how these methods are related and that whether ASEP is 
an exception or one can relate different one-dimensional 
driven-diffusive systems to the path models using these 
methods. On the other hand, one can also consider the 
driven-diffusive systems with periodic boundary condi- 
tions and investigate whether this connection still holds. 
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